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ABSTRACT 

Low-luminosity gamma-ray bursts (GRBs) with luminosity < lO'^^ergs"' probably consititute 
a distinct population from the classic high-luminosity GRBs. They are the most luminous 
objects detected so far within ~ 100 Mpc, the horizon distance of ultra-high energy cosmic 
rays (UHECRs), so they are considered to be candidate sources of UHECRs. It was recently 
argued that the energy production rate in UHECRs is much larger than that in gamma-ray 
photons of long GRBs measured by the Fermi satellite, which, if true, would challenge the 
view that GRBs can be the sources of UHECRs. We here suggest that many of the low- 
luminosity GRBs, due to their low luminosity, can not trigger the current GRB detectors and 
hence their contribution to the local gamma-ray energy production rate is missing. We find 
that the real local energy production rate by low-luminosity GRBs, taking into account the 
missing part, which constitutes a dominant fraction of the total amount, could be sufficient to 
account for the flux of UHECRs. Due to the low-luminosity, only intermediate-mass or heavy 
nuclei can be accelerated to ~ 10^" eV. We discuss the acceleration and survival of these UHE 
nuclei in low-luminosity GRBs, especially in those missing low-luminosity GRBs. At last, 
the accompanying diffuse neutrino flux from the whole low-luminosity GRB population is 
calculated. 
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Although the transition energy from galactic to extragalactic origin 
in high-energy cosmic ray spectrum remains inconclusive, there 
is a general consensus that cosmic rays with energy above the 
Greisen-Zatsepin-Kuzmin (GZK) energy of £ > £gzk = 4 x 10" 
eV, are of extragalactic origins. These UHE particles, whether 
they are protons or heavy nuclei, are attenuated when they are 
propagating in the intergalactic space. UHE protons with energies 
above the GZK cutoff undergo photopion interactions with cosmic 
microwave bacfcground (CMB) photons with an attenuation iength 
of ~ 100 Mpc. Coincidentally, heavy nuclei with energy > 50EeV 
will suffer from a strong photo-disintegration attenuation due to 
interactions with CMB and cosmic infrared background (CIB), 
with an attenuation length < 100 Mpc. As a result, sources 
producing UHECRs above E > Eqzk must be within ~ 100 
Mpc. Within this so-called GZK horizon distance, there are 
few sources that are powerful enough to be able to accelerate 
particles to energies ~ 10^" eV. The candidates include local 
active galactic nucleus (AGN) jets (e.g. iBiermarm & Strittmatteil 
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19871: iBerezinskv Gazizov. & Grigorievj |20 06|), local GRBs 

Wick, Derm er. & AtovarJ 



(Waxman 1995, 2 0041; IVietril . . , 

12004: iDermer & Atovanl l2006l: iMurase & Nagatakil 1200^ 
iMurase et alJ l2006l), and semi-relativistic hypemovae rem- 
nants jWang et alj bOOTa"). To accelerate UHECRs, the lu- 
minosity of the accelerators must satisfy the requirement 
L > 1.5 X lO'*^ergs-\r^/j3)(E/lQ^°eVf(Z/26)--, where F and /3 
are the bulk Lorentz factor and the velocity of the shock respec- 
tively, Z i s the nuclear charge number o f accelerated particles 
dWaxmaiJ l2005l : |f arrar & Gruzinovl 1 20091 : iLemoine 8i Waxmml 
I2OO9I) . The brightest sources within this distance are nearby 
low-luminosity GRBs (LLGRBs) associated with hypemovae, e.g. 
GRB980425 with peak luminosity L ~ 5 x 10'*''ergs"' at a dis- 
tance 40 Mpc (associated with SN1998bw) and GRB060218 with 
L ~ 5x lO^^ergs"' at distance 140 Mpc (associated with SN2006aj). 
They are much dimmer than their high-luminosity brothers, which 
are, however, detected at high redshifts. As a special subclass 
of GRBs, these nearby GRBs/hypern ovae hav e been proposed 
to be candidate sources of UH ECRs jMurase & Nagataki 200^ 
IWang et alJl2007ah . lWang et alj feOOvQ) suggest that GRB060218- 
like LLGRBs may arise from the breakout of semi-relativistic 
(r ~ 2), radiation-dominated shocks from the progenitor stars of 
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hypemovae or dense stellar wind surrounding the progenitor stars. 
In this scenario, UHECRs can be accelerated at the forward shock 
formed when the semi-relativistic eje cta is expanding in th e stella r 
wind. In the scenario proposed by iMurase & Nagatakil ( |2006|) , 
LLGRBs are thought to arise from internal shocks in relativistic 
outflow with r ~ 10, similar to the case of high-luminosity GRBs, 
and UHECRs are accelerated by the same internal shocks. 

GRBs are discovered preferentially at high redshifts (z ^ 0.5) 
with an isotropic-equivalent energy of Ey > lO^'^'^^erg, released in 
a few tenth of seconds to a few tens of seconds. So far only sev- 
eral nearby GRBs within ~ 200 Mpc are detected. GRB 980425, 
associated with hypemova 1998bw, is the first-found peculiar burst 
detected by BeppoSAX and BATSE at a distance of 38 Mpc, with 
an isotropic-equivalent total emitted energy of only ~ 9.3 x 10"*^ erg 
in 1-lOOOOkeV bancQ and a duration of ~35s dGalarna et al.ll 19981 : 



lKulkarmeta"i]|l998 



Plan etal.|[l999l : iKaneko et all 120071) . GRB 



031203, associated with hypemova 20031w, was detected by INTE- 
GRAL at z = 0.105 with an isotropic-equivalent bolometric energy 
of ~ 1.7 X 10^°erg and a duration of ~37 s jMalesani et aT]|2004i : 
IProchaska et"ai]|2004l : [Kaneko et alj|2007l) . GRB 060218, associ- 
ated with hypemova 2006aj, was detected by Swift at a distance of 



140 Mpc, with an isotropic-equivalent energy of ~ 4.3 x 10"" 



erg 



and a duration of ~ 2100s (Campana et al. 20061 : iMirabal et al.l 
l2006l : IPian e?a l. 2006; Kaneko et alj |2007|) . GRB100316D, as- 
sociated with SN 2010bh, was detected by Swift at distance of 
260 Mpc, with an isotropic-equivalent bolometri c energy larger 
than 5.9 x 10"*^ erg and a duration about 1300s ( Chornock et al.l 
l20ld : IStarling et alj|20in. All these nearby GRBs have low lu- 
minosity and are therefore named low-luminosity GRBs (LL- 
GRBs). Due to their proximity, the inferred intrinsic rat e of LL- 
GRB s is, however, much higher ( Soderberg et al. 2006; Li ang et al] 
l2007l : iGuetta & Delia Vallell2007l : IDai 2009) than high-luminosity 
GRBs (HLGRBs). By checking whether the event rate of LL- 
GRBs is consistent with a natural extrapolation of HL GRBs to 
low luminosity in a coherent luminosity function (LP), it is foimd 
that LLGRBs l ikely form an intrinsically distinct population from 
HLGRBs (e.g. iLiang et al."2007': 'C uetta & Delia ValldlioOTl : lOail 
l2009l : IVirgili, Liang. & Zhang,2009) . 

Recently, Eichler et alj ( 120101) find that the energy produc- 
tion rate in UHECRs above 4EeV is about 10^^ times larger than 
that contained in gamma-rays recorded by Fermi from long GRBs. 
This appears to be an overestimate, as later calculation bv WaxmanI 
( I2OIQ) shown that the energy production rates in gamma-ray pho- 
tons of GRBs and extragalactic UHECRs are comparable if some 
factors are taken into account carefully. The discrepancy may be 
caused by a combination of three factors: 1) Eichler et al. (2010) as- 
sumed that the extra-galactic CR production rate is about 10 times 
larger than the UHECR production rate, arguing that the genera- 
tion spectrum of the extra-galactic cosmic rays extends to energy 
smaller than 10"* eV (note the spectrum is dn/ds oc £"^'), while 
the estimate by Waxman (2010) only includes UHECRs above 10" 
eV; 2)they assumed different transition energy at which the flux of 
extra-galactic cosmic rays dominates over the flux of the Galactic 
cosmic rays. Eichler et al. (2010) assumed that extra-galactic cos- 
mic ray flux dominates over the Galactic component at 4 x 10'** eV, 
while Waxman (2010) assumed that it is dominated above 10" eV; 



' hereafter, unless otherwise specified , the "energy" or "luminosity" refers 
to this isotropic-equivalent bolometric (in 1-lOOOOkeV band) one 



3) Waxman (2010) also pointed out that the gamma-ray production 
rate in Eichler et al. (2010) is under-estimated, arguing that many 
distant/faint bursts are missed by Fermi/GBM. 

Motivated by this problem, we revisit the scenario of LLGRBs 
as the accelerators of UHECRs and study whether the energy bud- 
get in LLGRBs is sufficient. Because of the low luminosity of LL- 
GRBs, a lot of them may not trigger the detector and have therefore 
eluded detection. However, these dim GRBs may be still powerful 
enough to accelerate UHECRs. 

For these LLGRBs, only inter-mediate mass or heavy nuclei 
can be accelerated to energies ~ 10^" eV. The composition of the 
observed ultra high energy cosmic rays remains disputed. Recent 
observations of the Pierre Auger Observatory (PAO) show a tran- 
sition in the maximum shower elongations < X^a ^ > and in their 
fluctuations RMS(Xn,ax) between 5EeV and lOEeV j Abraham et al.l 
[2OIO). These transitions are interpreted as reflecting a transition 
in the composition of UHECRs in this energy range from protons 
to intermediate mass nuclei. However, this claim depends on the 
poorly-understood hadronic interaction models at such high ener- 
gies. There is a long-lasting tension between the spatial UHECR- 
Active Galactic Nuclei (AGNs) correlation suggesting protons, and 
Auger results of X^^^ and < > suggesting a significant heavy 
element component for UHECR in the same energy range. Recent 
update results of PAO show that the significance of the correla- 
tion between UH ECRs and nearb y extragalactic matter is decreased 
jThe Pierre AUGE R Collabora tion et alj|20ldl) , which relives this 
tension to some extent and allow the possibility of a heavier com- 
position. 

To show whether LLGRBs can be sources of UHE nu- 
clei, one also needs to know the survival probability of UHE 
intermediate-mass or heavy nuclei in the sources and the ori- 
gin of these nuclei. For the semi -r elativi stic hypemova sce- 
nario, 'W ang. Razzague. & MeszarosI j200 8*) have shown that 
intermediate-mass or heavy nuclei can easily survive in the hyper- 
nova remnant shocks as the shock size is typically very large. In this 
scenario, one would also expect a natural origin of intermediate- 
mass to heavy nuclei since hypernova remnants are expanding in 
the stellar wind of the progenitor star of hypemovae, in which inter- 
mediate mass n uclei, such as O and C , are enriched. For the internal 
shock scenario, IMurase et alj ( I2OO8I) discussed the survival proba- 
bility of heavy nuclei in LLGRBs for a certain set of parameter val- 
ues of the luminosity, shock bulk Lorentz factor and peak energy 
of the photon spectrum. However, these parameters may vary for 
different LLGRBs and there are some inherent correlations among 
them. We will study their effects on the survival of UHECR nuclei, 
especially in those dim LLGRBs that do not trigger the detector. If 
these non-triggered LLGRBs are also able to inject UHECRs into 
the universe, they will contribute to the energy production rate in 
UHECRs, but not to that in gamma-ray photons recorded by detec- 
tors. 

The paper is organized as follows. In §2, we estimate the ratio 
of gamma-ray energy production rates between the missing local 
LLGRBs and the observable local LLGRBs, as well as the local 
gamma-ray and UHECRs energy production rates for two lumi- 
nosity functions. We discuss the acceleration and survival process 
of UHE heavy nuclei in sources, especially in those LLGRBs that 
may be missed by the detector in §3. In §4, we calculate the ac- 
companying neutrino flux, contributed by all LLGRBs in the whole 
luminosity range. We give our conclusion in §5. Throughout the 
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paper, unless otherwise specified, we use eV as the unit of particle 
energy and use c.g.s units for other quantities and denote by 2, the 
value of the quantity Q in units of 10"^. 



2 THE LOCAL ENERGY PRODUCTION RATE BY 
LLGRBS 

To know the local energy production rate by LLGRBs, one needs 
to know the luminosity function (LF) of LLGRBs and then sum 
the contributions by all LLGRBs over the whole range of luminos- 
ity. Considering that LLGRBs probably form a distinct population 
from HLGRB s, here we a dopt the LFs of the broken power-law 



1000s. We 



form gi ven by Liang et ilOOi) andlDail ( |2009|) . The first one has 



a form jLiang 



Liang et aij 
et altlOGTh 



dN 



L 



(1) 



where po is the local event rate of LLGRBs inferred from the 
observed ones, and Oo is a normalization constant to guarantee 
the integral over the luminosity function being equal to the lo- 
cal event rate po. In this luminosity function, the total local LL- 
GRB rate is insensitive to the minimum luminosity, but fixed by 
the break lum inosity Li,. We tak e the best fit values for these pa- 
rameters from lLiang et al.l ( l2007t) in our following calculation, i.e. 
Po = 325Gpc"'yr"', Lb = iC^^ergs"', ai = 0, and = 3.5. 
Another LF for LLGRBs is suggested as (e.g. lDail2009h 



dN 

11 



--Po 



(2) 



which describes HL and LLGRBs together in one form with dif- 
ferent power-law index, joined at the break energy. Here we take 
Po = 1.7Gpc-3yr-', Lb = 5 X 10'*'*ergs-', ai = 2.3 and ci-j = 1-27, 
which are suggested bv iDail 12009) as the best fit values. In this 
LF, the local event rate of LLGRBs is sensitive to the minimum 
luminosity. The total event rate of LLGRBs is hard to know since 
the majority of LLGRBs at the low-luminosity end could be missed 
from detection. In principle, the total event rate of LLGRBs must be 
lowe r than the total rate of type Ib/c supemovae, ~ 2x 10'*Gpc ~^yr~' 
(e.g. ICappellaro. Evans, & Turattdfl^ : lOahlen et alj|2004h . The 
fraction of Ib/c supernovae (SNe) which have relativistic outflows 
is a somewhat uncertain number. Radio observations of a large sam- 
ple of Ib/c SNe suggest that less than 10% of Ib/c SNe are associ- 
ated with GRBs ( Soderberg et al, 2006). We assume two represen- 
tative values for the ratio of the local event rate of LLGRBs to that 
of Type Ib/c SNe, i.e. ^ = 1% and <f = 10% in the following calcu- 
lation. With these ratios, one can get the minimum luminosity L„,j„ 
for the LF, as shown in Table 1 . 

With the LFs given above, we can now calculate the local 
gamma-ray energy production rates by LLGRBs 



Wy(0) 



(3) 



where Ey{L) is the isotropic equivalent gamma-ray energy out- 
put for LLGRBs with luminosity L. The relation between Ey and 
L is unknown for LLGRBs due to the small sample detected 
so far. We assume Ey oc L*, where we take a wide range for 
k, i.e. k e (0, 1). We take the isotropic gamma-ray energy in 
GRB 060218 and GRB 1003 16D as a reference value, i.e. Ey = 



10"" ergs-'. For the LF 



LTgo = 10 erg, where L = 10 ergs" and Tgo 
se t the upper l imit o f integral at Ln,a). 

in iLiang et al] ( |2007|) (hereafter, LFL), there is no constraint on 
L,ni„ since the total local LLGRB rate is insensitive to the min- 
imum luminosity. Setting Lmin to 5 x 10*^ ergs"', which is suffi- 
ciently low to contain all the LLGRBs that have been observed 
to date, we have Wy(0) - 3.25 x lO^'ergMpc^Vr"' for '^=0 and 
2.72 X 10'''ergMpc"^^yr"' for k=l. The results are comparable for 
the case of L^]„ = 5 x lO ^'^ergs"', as shown in Table. 1. On the 
other hand, for the LF of iDail l l2009i) (hereafter, LFD), depend- 
ing on the minimum luminosity Ln,i„ of the luminosity function, 
Wy(0) is in the range of (2 - 20) x lO^^ergMpc-^r"' for k^O and 
(8.42 - 14.9) X lO^'ergMpc-^r"' for k = 1. The results of Wy{0) 
are given in Table 1 . These values are larger than that obtained in 
lEichleretai] ( l2010h by one to two orders of magnitude. Further- 
more, the fact that the energy production rate in the k = I case is 
less than that in the k = case in the lower L^^y, case indicates that 
most gamma-ray energy production rate is contributed by LLGRBs 
with relatively low luminosities (e.g. L < 10'*'' ergs"'), which may 
be difficult to be detected by current GRB detectors. 

We now discuss how much the missing LLGRBs contribute to 
the local energy production rate. To be detected by Fermi/GBM, 
a GRB located at a distance D requires to be brighter than a 
limit luminosity Lnn, = 4;rD^5GBM/S(«H «2)> where 5gbm = 
1.75 X 10"^ergcm"^s"' is th e sensitivity of Fermi/GBM in its 
burst trigger band ( 50-300keV, lvon Kienlin et al.l2004l ; lBandl2008l: 
llmeritoetal] |2008') and B{ei,e2) is the energy fracti on in the de- 
tector frequncy window. B(ei, 62) is given by CBloom. Frail, & S^H 
l200lHlmerito et al.ll2008h 



f'' Bn(E)ds 
j^^ En(E)dE 



(4) 



where e\ and 62 are the upper and lower threshold energy of the 
detector, while £1 and E2 are upper and lower limit for the bolo- 
metric gamma-ray spectrum. For Fermi/GBM, ei = 300keV and 
£2 = 50keV. To be consistent with the inferred bolometric lumi- 
nosities of those observed LLGRBs, as mentioned in Sec.l, we 
set El = lOOOOkeV and & = IkeV. n{E) is the LLGRB prompt 
photon spectrum, where e is the photon energy. We assume a bro- 
ken power-law spectrum similar to that of HLGRBs, expressed by 



dnIdE = nt,(E/Byi,) 



for e < Syb and dn/ds - n^,(s/Ey^,) 



for 



E > e^b- We take /?i ^ 1 at energies below the break and /?2 - 2 
above the break. Then the gamma-ray energy production rate by 
those LLGRBs that are missed by Fermi/GBM can be estimated by 



Wy. 



1 



r-Omux /-i|,„, 



4nD-Ey{L)^dLdD, (5) 



where Li,„, = 2.1x 10'"'(D/100Mpc)VS(ei, e2)ergs-' . Since 40Mpc 
is the distance of GRB 980425, the closest LLGRB ever detected, 
while 200Mpc is approximately the attenuation length for iron nu- 
clei of lOOEeV, we take £)„,]„ = 40Mpc and D^.^ = 200Mpc as the 
lower limit and the upper limit of the integration respectively. Due 
to that Ey], varies with luminosity, S(ei, £2) is a function of the lumi- 
nosity. We assume that the Yonetoku relation holds for LLGRBs as 
in the case of GRB 060218 (Yonetoku et al. .2004) . i.e. e^b 
and take £^b ■ 



lOL^^'keV. We define 

H = Wy.nJ(Wy - Wy.n.), 



(6) 



as the ratio between the gamma-ray energy production rates in 
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Liang et al, 2007, L,,,i„=5x 1 0^ergs" 

Liang et al. 2007, L^i„=5x1 o^^ergs"^ 

- - Dai 2009, e=10%, L„i„= 1 .79x 1 0"ergs"' - 
Dai 2009, J=l%, L„;„=l,09xlO*'ergs" 
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Figure 1. The ratio of the missing to the observable gamma-ray production 
rate % changes with the assumed power-law index of the Ey - L relation 
k {Ey = 10'''L^yerg). Except in the f = 1% case for LFD, the missing 
gamma-ray energy production rate is several times larger than the observ- 
able one. 

missing local LLGRBs and that in the observable ones, and present 
its values in different cases in Table. 1. We find that the gamma- 
ray energy production rate by missing local LLGRBs constitutes a 
dominant fraction of the total rate. In some cases, it is larger than 
that produced by the observable LLGRBs by one order of magni- 
tude. Note that in the ^ = 1% case of LFD, that K becomes small is 
due to an unrealistic high Lmi„, which is higher than the luminosity 
of detected LLGRBs (e.g. GRB 980425, GRB 060218). We show 
how K changes with k in Fig. 1. 

Although the obtained gamma-ray energy production rate is 
highly dependent on the typical value of Ey used (i.e. we used 
Ey = 10^''erg for a LLGRB with luminosity of 10"*' ergs"'), which 
could have large uncertainty, the ratio K is only dependent on 
B(ei,e2), or to be more accurate, the assumed photon spectrum. 
We also check some other common combinations of /JnySa and e^j, 
and present the results in Fig. 2. In the region that L < 10'*'' ergs"', 
where the missing LLGRBs reside, the value S(ei, 62) is generally 
smaller than 0.3. It is mainly because that the energy window of 
the detector is so narrow that the spectral peaks for LLGRBs are 
typically outsi de of the windo w. 

Recently, lEichled l [20Tll) have shown that the contribution of 
dim or undetected GRBs to the total all-sky GRB energy flux is 
small, based on that the cumulative fluence plateaus at a fluence 
of ~ 10"'ergcm"-, which is much higher than the minimum flu- 
ence for both BATSE and Swift. However, we note that the trigger 
threshold for SwiftfBAT and Eermi/GBM is about 10"*ergcm"^s"' 
and LLGRBs typical have long duration of T ~ lO^s, so LLGRBs 
with fluence lower than ~ 10"'ergcm"^ would not trigger the detec- 
tor. Therefore, the plateau only reflects that those GRBs at the low 
luminosity end of the HLGRB population contribute little to the 
cumulative fluence (which is actually consistent with the LFs of 
HLGRBs), irrespective of the LLGRB population. As a distinctly 
separate class of bursts at very low flux levels, LLGRBs could make 
a significant contribute to the cumulative fluence only if the detec- 
tor is significantly more sensitive than Swift/BAT and Fermi/GBM. 

In the internal shock scenario for the gamma-ray emission of 
LLGRBs, UHECRs are accelerated by the same shocks. So one can 



Figure 2. Relation between bandpass ratio B(e\,e2) and luminosity with 
different combinations of parameters of photon spectrum. See text for more 
discussion. 

estimate the energy production rate in UHECRs, assuming some 
equipartition factors for electrons and protons (nuclei). Assuming 
a flat spectrum dh/ds oc for UHECRs, we can now estimate 
the local UHECRs energy production rate per logarithmic energy 
interval by LLGRBs 

2dh 1 . 
Wuhecr(O) = e -j-luHECR = 7— ; tIVcrCO) 

de in(eA,max/£A,min) 



ln(eA,max/£A.min) Q 

where Wcr{0] is the local energy production rate for cosmic ray 
over the whole energy range, £A,min and SA.max are, respectively, the 
minimum and maximum energies of cosmic rays. Though the value 
of e^.rain is unknown, it is usually assumed that ln(eA.n,a,x/£A,miii) - 
10. Here and 64 are, respectively, the equipartition factors for 
electrons and protons (nuclei). We adopt the usual values, = 
Eg = 0.1 and Ea = 1 - - (^b is the equipartition factor 
for the magnetic field), in the following calculation. The results 
of fi^^luHECR are given in Table 1 for the two LFs considered. 
These values are not far from the measured value of the local en- 
ergy production rate of UHECRs, (5 ± 2) x 10'*^ er g Mpc'^yr ' (e.g. 
iKatz. Budnik. & Waxmaiill2009l : IWaxmanll2O10l) . so we conclude 
that local LLGRBs remain viable sources for UHECRs in terms of 
the energy production rate. 

In the semi-relativistic hypernova scenario for the gamma-ray 
emission of LLGRBs, UHECR acceleration occurs in the hyper- 
nova remnant shock region, which is different from the gamma- 
ray production region. Nevertheless, the energy production rate in 
gamma-rays is still a useful indicator, as it reflects a minimum en- 
ergy production rate in the kinetic energy of the semi-relativistic 
ejecta, from which the energy of UHECRs is tapped ultimately. 



3 UHE NUCLEI FROM LLGRBS 
3.1 Acceleration of UHE nuclei 

Now we consider the acceleration of nuclei by semi-relativistic hy- 
pernova remnant shocks or internal shocks in relativistic shells. For 
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Table 1. The ratio of gamma-ray energy production rates between the missing local LLGRBs and the observable ones, estimated gamma-ray energy production 
rate of local LLGRBs and the inferred UHECRs energy production rate in per logarithmic energy interval. See text for more discussion 



Luminosity Function L„,in(10 ergs ) 



■R" Wj,(0)'"(10''3ergMpc-3yr"') IVuHECRCOfClO^'ergMpc-^yr-') 

k=0 k=l k=0 k=l k=0 k=l 



Liang et al. (2007) 0.05 7.1 2.2 3.25 2.72 2.60 2.18 

Liang et al. (2007) 0.5 4.0 1.8 3.25 3.95 2.60 3.16 

Dai (2009) (po = 2000Gpc-3yr-i) 0.18 12.4 1.1 20.0 14.9 16.0 11.9 

Dai (2009) (po = 200Gpc-3yr"') 109'' 0.89 0.21 2.00 8.42 1.60 6.74 



" : The ratio of gamma-ray energy production rates between the missing local LLGRBs and the observable ones. 
: The total (the missing + the observable) gamma-ray energy production rate of local LLGRBs. 
: The inferred local UHECRs energy production rate in per logarithmic energy interval, see Eq. (7) 

: Since the luminosities of GRB 980425 and GRB 0602 1 8 are less than this value by a factor of 2—4, this case may be largely dilferent from the reality. We 
just show them here for reference. 



Fermi shock acceleration, the characteristic acceleration timescale 
is face = '^^IPlh' where A = t}R[^ is the scattering mean free path, 
Rl is the Larmor radius of particles and /3shC is the shock speed. 
Here ij ~ a. few, describes the ratio between the acceleration time 
and Larmor time (77 = 1 corresponds to the efficient Bohm diffu- 
sive shock acceleration). For a nucleus of mass Am,, and charge Ze, 
the time needed to accelerate it to the observed energy of ea in the 
comoving frame is t.^^^ = jjeaI {ZeBY0l^c). Cooling of particles re- 
stricts the maximum energy. Adiabatic expansion and synchrotron 
emission are two important cooling mechanisms. The cooling time 
of the former one is given by 4d - Uyn - Rsh/^PihC, where i?sh is 
the radius of shock acceleration site and F is the bulk Lorentz fac- 
tor of the shock, while the latter one is t.y„ = — ; — 52il£ . vvhere 

(Tr is Thomson scattering cross section for electrons and tr7-(^)Z'* 

"'a 

is the corresponding cross section for nucleus with Ze charge and 
Anip mass. 

3.1.1 Semi-relativistic hypernova scenario 

Whether particles can be accelerated to ultra-high en- 
ergies in the hypernova s c enario has been discussed in 
IWang. Razzaque. & MeszarosI ( |2008[) . in which a distribution 
of the ejecta energy with ejecta velocity has been assumed. 
Here for simplicity we only focus on the semi-relativistic 
part of the ejecta, whose velocity is Tp > 0.5. According to 
the simulation of hypernova explosio n such as SN1998bw by 
iMacFadven. Wooslev. & HegeJ bOOll) , the isotropic-equivalent 
kinetic energy of ejecta is roughly constant at angles larger than 
20°, so we assume a spherical hypernova ejecta expanding into 
the circum-stellar wind medium. Particles are accelerated in 
the region where the ejecta is freely expanding before being 
decelerated by the swept-up circum-stellar medium. The size of 
this free-expansion phase region for ejecta of a particular velocity 
jSshC and kinetic energy is 

Ahn - 4 X 10"£,,M (FyS,h)"'M!f v-,;/-cm, (8) 



which is much larger than the radius at which gamma-ray pho- 
tons of LLGRBs are produced (Waiig et al- 2(X)7bJ, where M = 
iO^^M-sMgyr"' is the wind mass loss rate, whose average value 
is 3 X lO'^Mo vr'' for WR stars, and v„, = lO^ Vw.3 kms"' is the 
wind velocity Jwinislll99ll : IChevaUer & Lilll999h . During the free 
expansion phase, the magnetic field energy density is B^/Sn = 
2eBP„ (i?HE)c^/3si,, where eg = 0. 16b__i is the fraction of the equipar- 
tition value of the magnetic field energy and p„. is the mass density 
of the stellar wind at radius /?hn- The magnetic field at the free- 
expansion radius Run is 

B = 0.54^:,i?HN,,7ySshMi//v;'/^G. (9) 

From ?a;;c = 'dyn, the maximum energy is 

£A,max - ZeBR^^^Jri = 5.2 x lO^^-' (^) 4!-i/8sh'^'^';'f eV. 

(10) 

Note that the synchrotron loss of UHE nuclei in the semi-relativistic 
hypernova scenario is much lower than the adiabatic loss, so it is 
not considered here. 

3.1.2 Internal shock scenario 

In the internal shock scenario, the magnetic field B can be estimated 
by B- /Sn = egU = -^-Uy as all the energies of electrons are almost 
lost into radiation, where U is the comoving thermal energy density 
and Uy = L/4nR^^T^c is the comoving gamma-ray energy density 
for an internal shock radius of Rsh. So, 

^ ^ (^T fk = ' io'v^«:-.4f rr^*2-'G (11) 

where R^t, = 2T^c6t has been used and St is the variability time. 

The maximum nucleus energy is determined by equating /acc 
with the smaller one of the two cooling timescales fad and fsyn. We 
have 

£A,n,ax = 2.0 X 10^'77-'(:^)e;;lf 6^^!,L^/-r['/3„„eV (12) 
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for = fjy„, and 

(13) 

for fjicc = 4yn- From Eq. il2\ . one can derive a minimum luminosity 
for LLGRBs that can accelerate nuclei to ultra high energies, i.e. 

Z 



2.5 X 10^^el,ooEev'7'(^)"f.,-if^:-ir?/?rJergs-'. 



(14) 



According to Eq. ( I13K there is also a maximum luminosity for LL- 
GRBs that can accelerate nuclei to ultra high energies, i.e. 

'56' 46' 

(15) 

The steep dependence on the bulk Lorentz factors T results from 
the steep dependence of the magnetic field energy on T in the in- 
ternal shock assumption. If the luminosity is too high, the radiative 
cooling is so rapid that these nuclei can not reach ultrahigh ener- 
gies. As the minimum luminosity adopted in Sec. 2 is higher than 
^min.acc, We couclude that even those missing LLGRBs are powerful 
enough to accelerate UHECRs. Since the internal shock is mildly 
relativistic with /Jj,,, a 1, hereafter, we drop the dependence on Pi„, 
in our analytical calculation. 



where ctgdr = 1.45 x lO^^^Acm^ is the peak cross section of the 
Giant Dipole Resonance (GDR), p ^ 4 is the compression ratio of 
the hypernova remnant shock and ehe - leV is the characteristic 
energy of hypernova thermal photons. 



3.2.2 Internal shock scenario 

1. Photon Spectrum 

The photon spectrum of the prompt emissions of GRBs can be 
approximately described by a broken power-law. In the comoving 
frame of the relativistic outflow, it can be described as 



n{Sy) = ni, 



E < S) 

e' 



(18) 



where is the break energy in the comoving frame. As we did in 
Sec. 2, we take ySi = 1 and = 2 as typical values. At low energy 
region, in the framework of the internal shock model, there may be 
another spectral break, i.e. the synchrotron self-absorption (SSA) 
break. 

Now we estimate the SSA break energy for LLGRBs. The two 
characteristic break frequencies in the synchrotron spectrum in the 
comoving frame are given by 



3.2 Survival of UHE nuclei 

UHE nuclei will suffer from photo-disintegration and photopion 
loss by interactions with low-energy photons in the sources. 
For UHE nuclei, photo-disintegration is usually more important 
jWang. Razzague. & Meszarosll2008l : lAllard et alj|2008l) . Now we 
study whether UHE nuclei accelerated in LLGRBs can survive in 
the semi-relativistic hypernova scenario and internal shock sce- 



3.2.1 Semi-relativistic hypernova scenario 

There are two photon sources which could cause photo- 
disintegration of heavy nuclei: one is provided by hypernova ther- 
mal photons from radioactive elements of the hypernova ejecta, and 
another is provided by the synchrotron photons from the hypernova 
remnant shock. The free-expansion time for the semi-relativistic 
ejecta is 



™ 1.3x 10'£*.5i(rAh)"'M 



5 '^i,.',3 ^• 



(16) 



At earlier times when the fast semi-relativistic ejecta is decelerated, 
the flux from the hypernova thermal photons is expected to be dom- 
inated. We use the luminosity of SN1998bw as a representative for 
the hypernova luminosity of thermal photons. At time / ~ 100 days 
after the burst, the optical lumi nosity of SNI998bw d rops to the 
level of about Lhe ~ 10'"ergs ' jSollerman et alj|2002h . A nucleus 
of energy E = 10^" eV interacts with target photons with a thresh- 
old energy Sth <; 0.01(A/56)£'2o'eV. A rough estimate of the optical 
depth of photo-disintegration of UHE nuclei (losing one nucleon) 
due to hypernova thermal photons is 



T < o-aoRi- 



-)(^) 



(17) 



1.02X 10'V(p)r7^^ftj'e;'_^,e 



1-^47 



(7u 



if Hz 



and 



2.37 X 10''r^(5f2e^''_-,egi',-L:,f-(l + YT^Hz, 



(19) 



(20) 



where f(p) = 6(p - 2)/(p - 1), yim is the Lorentz factor of the in- 
ternal shock, and Y is the Compton parameter for inverse Compton 
loss. In the fast cooling regime (v^. < v^), Y = -i+ Vi+'*^<'/ a ~ q g 
for Eg = Eg = Since it is less than unity, we just neglect its influ- 
ence on the value of Vc- The SSA co efficient in fast cooling regime 
can be estimated by jWu et al.ll2003h 



-5/3 
3 



fMt) 



-(p+S)/2 



V < Vc 

Vc <V<V,„ 

V„, < V 



(21) 



where y,,, and jc is the minimum Lorentz factor and cooling Lorentz 
factor, «(, is the number density of electrons and cq ~ 15 is a con- 
stant. The SSA frequency is determined by kyAR^^ = 1, where AR^^ 
is the length of the shock region in the comoving frame. AR^o is re- 
lated to the column density of electrons E by 



(22) 



where N\ 



is the total number of electrons in the 



'shell ee(rinl-l) Tm^c- 

coUiding shell. Finally, we get the SSA frequency in the comoving 
frame, i.e. 



SSSA 



: o.oi2r7'^^/«&-"'"*Li5/^«eV 



(23) 



for Ec = 0.1, eg = 0.1, 7in( = 2 and p = 2.2. 

Since vssa > v„,, the SSA coefficient is oc e^^-^ for p = 2.2 
according to Eq. i2U . hence the spectral index below the SSA fre- 
quency is -2.6. Thus, the photon spectrum in the comoving frame 
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(essA/e"b)"'(e/£ssA)'',e < £ssa 



LuiTiinosit^^(erq/s) 



n(E) 



■CO \-l 
2 



ESSA <E<£ 



max ■ 



(24) 



2. Photodisintegration 



For an UHE nucleus with energy sa in tiie observer frame, tlie frac- 
tional photo-disintegration rate is 



-1 . ^ 1,^^, ^ 



I (i6crdis(e)e I lij^x ~n(x) (25) 



where = EA/TniAC^ = 1.9 x 10*(A/56) Ti 'E^aooEeV, 0-^,^(6) 
is the photo-disintegration total cross section, 6 is the photon en- 
ergy in the nucleus rest frame and ett, is the threshold energy. This 
rat e can be calculated numerica l ly wit h the cross section given 
in IPuget. Stecker. & BredekampI l ll976h . But as an estimate we 
can approximate the cross section mainly contributed by the giant 
dipole resonance (GDR) and show the numerical results later. The 
peak cross section due to GDR is iTgdr = 1.45 x lO^-'Acm^ and 
eoDR = 42.65A-''-2'MeV (for A > 4) with a width Agdr = 8MeV. 
We find that the energy of the photons that interact with UHE nuclei 
via GDR resonance is about e^^^ ~ 0.2ri(A/56)''™(£A.50Eev)"'eV 
in the comoving frame, which is larger than the SSA frequency, so 
the SSA process in the low-energy photons will have little influence 
on photo-disintegration process for these representative parameters. 
Approximating the integral by the contribution from the resonance, 
the fractional photo-disintegration rate is 



Uy cctgdrAgdr 



4e™ 



Aegdr 



(EA/EAbt-' 



E < EAb, 



(26) 



where Uy 



«b(ej 



^ )^[1 + In(e350kev/e7b)] ~ 4nb(e^t) is the en- 
ergy density of photons in the comoving frame in the energy 
window of SwiftfBAT and EAb = 0.5eGDR™AC^r^/£yb - 1-5 x 
lO'^AF^e^iJ iQijj.vSV is the nuclei break energy in the observer frame. 
For UHE nuclei, we have Ea » e^b^ then the elfective optical depth 
for photo-disintegration is 



Tdisi 



^dyn 
fdis.f 



: 0.043 



L47(A/56)°-" 

F|5f2£yb.l0ki;V 



(27) 



for jSi = 1. Thus we can see that a larger break energy, a larger 
bulk Lorentz factor or a smaller luminosity will be favourable for 
the survival of UHE heavy nuclei. We note that these quantities in 
GRBs may have inherent correlations. 

This motivates us to study in which kind of bursts UHE nuclei 
can survive, especially whether UHE nuclei can survive in those 
missing, dim LLGRBs, when these inherent correlations are taken 
into account. For the break energy Eyi,, we still assume that the Yo- 
netoku relation holds for LLGRBs, i.e. e^b oc L''^. As to F, assum- 
ing the break energy of the spectrum in the comoving frame are 
constant for all LLGRB^ we have F oc Eyi, and hence F oc L''^. 
The variability time scale 6t may depend on the central engine ac- 
tivity and we simply fix its value to 100s. Under these assumptions, 
we have t^isj «: L"^'^, and consequently, UHE heavy nuclei can 



Given that LLGRBs have lower break energies than HLGRBs, and on 
the other hand, LLGRBs are suggested to arise from less relativistic jets 
ISoderberg et"ai] j2006l) ; lToma et ajj ilOOik . the intrinsic range of the break 
energy in the comoving frame could be small. 



f-L, 50EeV, 
,-L. SOEeV, 
f-L, SOOEeV, 



"^aisj-EA. L=10 erg/s, E^a= 1 OkeV, r=10, (5t=100s^. 



10" 10' 10" 10" 

E. (eV) 



Figure 3. The fractional optical depth for photodisintegration of iron nuclei 
in the internal shock scenario. It is a function of the energy of nuclei (the 
solid line) and is a function of the luminosity of LLGRBs (dashed lines). 
The parameters used in the calculation are given in the figure for each lines. 
See text for more details. 



survive more easily in LLGRBs with relatively high luminosities. 
The main reason for this is that LLGRBs with lower luminosity 
may have smaller internal shock radii for a fixed 6t and hence a 
higher density of target photons. On the other hand, if we assume 
that F is a constant for LLGRBs, we will obtain r^i^s L^^^- In this 
case, UHE heavy nuclei can survive more easily in LLGRBs with 
relatively low luminosities. 

The SSA break may affect the photo-disintegration rate when 
^GDR ^ EssA> because the number of photons that can interact with 
the heavy nuclei via GDR is considerably reduced by the SSA pro- 
cess. This is applicable for nucleus at the highest energy (i.e. typi- 
cally above a few times lO^^eV). In the SSA regime, we have 



C6GDRO" GDR Agdr 
.2.84° AjIeI,^ 

The effective photo-disintegration optical depth is 

4"(A/56)'™ 



Td,s.f(£A) = 0.039 



1 J Ot^ fiyb.lOkeVtE/l.SOOEeV) 



(28) 



(29) 



If we assume e^b °^ L}^^ and F oc Sy^,, we will have r^is f oc L"''**, 
hence UHE heavy nuclei can survive from photo-disintegration 
more easily in LLGRBs with relatively low luminosities, if the en- 
ergy of the effective target photons falls below the SSA break en- 
ergy. 

In Fig. 3, we show the effective fractional photo-disintegration 
optical depth r^is f for iron nucleus, calculated according to Eq. l l25t 
and using the cross section described by the Lorentzian form in 
the energy range 6,1, < 6 < 30 MeV and a flat cross section for 
multi-nucleon loss in the energy ran ee 30MeV < e < 150MeV 
jPuget, Stecker. & Bredekamdll976lJKarakula & Tkaczvklll993l) . 
From the solid line, one can see that Tdis.f is almost independent of 
the energy of nucleus ea when ea > e^b (see the approximate ana- 
lytic estimate in Eq. l|27t). When the energy of the nucleus becomes 
larger than several times 10^° eV, the nucleus starts to interact with 
target photons with energies below the SSA break so that the effec- 
tive optical depth decreases with the energy of nucleus. Fig. 3 (the 
dotted, dashed and dash-dotted lines) also describes how rji, j of 
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Figure 4. The parameter space of luminosity L and Lorentz factor T in the 
internal shock scenario. The hatched area is the parameter space of L and T 
of LLGRBs in which iron nuclei of energy 50 EeV can both be accelerated 
and survive in the sources. The thick solid line represents the condition 
^dis.f = 1- The area between the two thick dash-dotted lines represent the 
parameter space in which iron nuclei can be accelerated to 50 EeV. Also 
shown is the Yonetoku's relation for GRBs (the thin dashed line). 



an UHE nucleus with a fixed energy depends on the luminosity of 
LLGRB. We can see that under the assumptions of £yb L''- and 
r oc £^1,, UHE nuclei tend to survive in LLGRB with relative high 
luminosity. The relation between r^is f and in the assumption that 
r is a constant for all LLGRBs is also presented in Fig. 3 for com- 
parison. In this case, UHE nuclei survive more easily in LLGRB 
with lower luminosity. 

In Fig. 4, we show the parameter space (denoted by the 
hatched area) of L and F in which iron nucleus can be acceler- 
ated to 50 EeV and can escape from the source with most of its 
initial nucleons preserved meanwhile. It shows that LLGRBs with 
larger L tend to satisfy both constraints in the internal shock sce- 
nario. Nevertheless, Tjisf ^ 1 occurs at L ^ 10'**ergs"', which is 
sufficiently low that most LLGRBs we are concerning are within 
the safety limit. 



4 NEUTRINO BACKGROUND FROM LLGRBS 

UHE nuclei will interact with low-energy photons in the sources 
and produce high-energy neutrinos. Hence neutrino detection is a 
useful tool for probing the acceleration of UHECRs in sources. In 
the semi-relativistic hypernova scenario, the production of neutri- 
nos is mainly through the photomeson in teractions o f UHE nuclei 
with hypernova thermal photons (Wang et alj|2007ai) . It has been 
found that the diffuse neutrino flux contributed by semi-relativistic 
hypernova is too low to b e detected by cubic kilometer detectors 
such as Icecube jWang et al. 2007a). 

In internal shock scenario, the production of neutrinos is 
mainly through the photomeson interactions of UHE nuclei with 
low-energy prompt photons in LLGRBs. The fractional energy loss 
rate of a nucleus with energy in the comoving frame due to pho- 



tomeson productions is 
- 1 1 

?mes(rA) = — I — lir 

JA at 



2ri Jf,,, J^i 



dxx ^n{x) 



(30) 



where e is the photon energy in the rest frame of the nucleus, 
6,h = 0.15GeV is the threshold photon energy for photomeson in- 
teraction, tTmcsCe) is the cross section of photomeson production 
and ^aCc) is the average fraction of energy lost to secondary pi- 
ous. Above the threshold energy, the main contribution of pion 
production is due to the A resonance at ea = 0.34GeV, for which 
the cross section ca n be approximately by a Lorentzian form (e.g. 
iMticke et al ] |200(]I) . Neglecting the nucle ar shadowing effect in 
photoproduction jMichalowski et al.lll977[). the peak cross section 
is Q-A - a-^A ^ 4.1 X lO'^^Acm^ jAnchordogui et alj l2008l: 



iMurase et aObOOSh . 



where apy is the peak cross section for pro- 
tons. As ^aC^) ~ ^p(f)/A, cT/i^A is independent of A. The energy 
of the photons that interact with UHE nuclei via A resonance is 
~ 3.8Fi(A/56)(ea,50Ecv) 'eV in the comoving frame . Since this 
energy is much larger than the SSA break energy, the SSA process 
has little effect on the photomeson process, even for the nucleus 
with energy of a few times lO'^eV. Approximating the integral by 
the contribution from the resonance, the total fraction of energy lost 
by nuclei to pions is 



^dyn 
^ mes 



Tct„ 



0.004 



£yb,10kevF'J(5/2 



(31) 



So f„(EA) depends on the inherent relations among F, Syt, and L. For 
Eyi, OC L''- and F oc Ey^,, we have f„(EA) oc L^^'^. But in the case that 
the bulk Lorentz factor is a constant with F = 10, fA.EA) °^ i-''^- For 
LLGRBs with luminosity in the range of L ~ 10**"'" ergs"' , f„ <ii I 
for both cases, which indicates that most heavy nuclei suffer from 
negligible energy loss due to photomeson production process. 

The diffuse neutrino flux contributed by all LLGRBs can be 
obtained through the integration over the whole luminosity range 
of LLGRBs, which is 



odNy, ^ 1 c 

SZ (Ev) = 

' dEy AAnHo 



L 



dEA 



dL 



(32) 

where /- ~ 3 is the correction factor for the contribution from 
high redshift sources and f„ is the suppression factor on the 
neutrino flux due to pion cooling i Rachen & MeszarosI 1 19981 : 



iRazzaque, Meszaros. & Waxmanll2004l ; IWang & Daill2009l) . Since 
the synchrotron cooling dominates the cooling process, we have 
& = T";'/(Cyn + '^^''>' where t^yn.,, = l2L'^^T]Stl{E„,,E^yy^s is the 
synchrotron cooling time and t„ = 2.6x lO^^y^ = 186£;rjEcvS is the 
lifetime of pions. The lower limit of the integral L^i„(EA) is given 
in Table. 1. The results obtained numerically are shown in Fig. 5 
for the two LFs. It shows that only in the case of high local rate 
of LLGRBs, there is a chance that diffuse neutrinos from LLGRBs 
could be detected by cubic kilometer detectors such as Icecube. 



5 DISCUSSIONS AND CONCLUSIONS 

The origin of heavy or intermediate-mass nuclei in LLGRBs re- 
mains to be an open issue. In the hypernova remnant scenario, 
cosmic ray particles originate from the material swept-up by the 
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Figure 5. The diffuse neutrino bacftground from LLGRBs. The top panel is 
for the LF in Liang et al. (2007), whereas the middle and the bottom ones 
are for the LF in Dai (2009), but for a local rate of po = 200Gpc"^yr"' and 
po = 2000Gpc"'yr"' respectively. The reference values of e^i,, F and iJy jso 
in the scalings are lOkeV, 10 and 10'''erg respectively for a LLGRB with lu- 
minosity of lO^'^ergs"' . We also plot the IceCube sensitivity to diffuse neu- 
trino fluxes in 4 years operation I Spiering 2008) and the Waxman-Bahcall 
bounds (Waxman & Bahcall 1999). 



semi-relativistic hypemova shock front. The progenitors of hyper- 
novae are thought to be Wolf-Rayet stars, as the spectral type of 
the discovered supernovae in these events is typically Ic. These 
stars are stripped of their hydrogen envelope and sometimes even 
the helium envelope, so the circum-stellar wind is rich of heavy 
elements, such as C and O. Thus, heavy or intermediate-mass 
UHE nuclei may naturally originate from the element-enriched 
stellar wind of Wolf-Rayet stars in the hypernova scenario. In WC 
type Wolf-Rayet stars, the C abundance is Xq ^ 20 % - 55% 
(by m ass) and the O abundance is Xq ^ 5% - 10% jCrowthe^ 
I2OO7I) . In WO type Wolf-Rayet stars, O and C abundances are 
even highe r, with Xc ^ 40% - 55% and Xq ^ 15% - 25% (e.g. 
Kingsburg h, Barlow. & Storevll995h . The abundances of these ele- 
ments are clearly much higher than the solar values. In the internal 
shock scenario for LLGRBs, the origin of nuclei is, however, quite 
uncertain. The uncertainty lies in whether there are nuclei entrained 
into the relativistic jet during the formation of the jet out of the col- 
lapsing core of massive stars. The temperature of the accretion disk 
resulted from the collapsing core needs to be sufficiently low that 
nuclei in the disk will not be photo-disintegrated by the disk ther- 
mal photons. This condition may be satisfied if the accretion rate 
for LLGRBs is sufficiently low. Another possibility for the origin 
of nuclei in the internal shock scenario is that nuclei present in the 
surrounding stellar envelope are entrained into the jet while jet is 
propagating through the stellar envelope. The validity of this possi- 
bility needs a detailed numerical simulation of the jet propagation 
in LLGRBs. 

In this paper, we have suggested that many dim local low lu- 
minosity GRBs may be missed by Fermi/GBM and these missing 
LLGRBs could make a dominant contribution to observed flux of 
UHECRs. If true, it can relieve partly the discrepancy between the 
energy production rate in UHECRs a nd that in gamma-ray photons 
recorded by Fermi/GBM, as raised by Eichler et al.l ( l20I(]|) . We first 
calculate the energy production rate in gamma-ray photons by LL- 
GRBs as a separate population from high-luminosity GRBs and 
estimate how large those missing LLGRBs contribute to the total 
gamma-ray energy production. In the calculation, we take two dif- 
ferent LFs for LLGRBs that were proposed in the literatures. We 
find that the gamma- ray energy production rate by LLGRBs for 
both LFs are one to two ord ers of magnitude larger than that esti- 
mated bv lEichler et alj JioTOi) . The missing part of the energy pro- 
duction rate in gamma-rays could account for a fair proportion or 
even the vast majority of the total one and hence could be much 
larger than the observable part. It should be noted that our results 
depend on many assumptions, such as the form of LF, the spectrum 
of LLGRBs and some inherent relations among parameters of LL- 
GRBs (e.g. Ey-L relation, Ej^ - L relation). To get a more accurate 
results, further studies in these relations with larger statistics are 
needed. 

There are two scenarios for the UHECR acceleration in LL- 
GRBs, one is the semi-relativistic hypemova scenario where UHE- 
CRs are accelerated by the semi-relativistic hypernova shock ex- 
panding into the circum-stellar wind and another is the internal 
shock scenario where UHECRs are accelerated by internal shocks 
within the variable relativistic outflow. We find that, in both scenar- 
ios, only intermediate-mass or heavy nuclei could be accelerated 
to energies above 10^" eV, which is consistent with the recent find- 
ings of heavy UHECR composition by R^O. We have also studied 
whether heavy nuclei could survive from photo-disintegration by 
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low-energy photons in the sources. In the semi-relativistic hyper- 
nova scenario, UHE nuclei can survive easily as the large size of 
the acceleration region leads to a low density of target photons. On 
the other hand, the survival probability of UHE nuclei in internal 
shock scenario depends on the inherent relations among the quan- 
tities of LLGRBs, such as e-yi,, T and L. In spite of the assumptions 
in these relations, we find that in the luminosity range of LLGRBs 
that we are concerning, UHE heavy nuclei can retain most of their 
nucleons before escaping from the sources. Finally, the accompa- 
nying neutrino flux is calculated. The energy loss through photo- 
pion production is negligible for UHE nuclei and hence the difiiise 
neutrino flux from LLGRBs can hardly be detected by current de- 
tectors, except in the case of a high local rate of LLGRBs. 

Our calculations indicate that local LLGRBs remain viable 
sources of UHECRs in terms of the energy production rate. Fu- 
ture GRB detectors, such as EXISTfHEl, could cover a wider 
burst trigger band (5-600keV) with an advanced sensitivity of 
1.5 X lO-'ergcm-^s-' (Imerit o et aL 2008.) . £X/Sr/HET will in- 
crease the LLGRB detection rate by a factor of tens compared to 
SwiftfBKI and Ferm;/GBM l ilmerito et alj|2008l) , thus much more 
LLGRBs could be detected. It would help us to learn more about 
the nature of LLGRBs and examine more carefully the connection 
between LLGRBs and UHECRs in future. 
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